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Abstract 

Non-perturbative studies of the thermodynamics of strongly interact- 
ing elementary particles within the context of lattice regularized QCD are 
being reviewed. After a short introduction into thermal QCD on the lat- 
tice we report on the present status of investigations of bulk properties. 
In particular, we discuss the present knowledge of the phase diagram in- 
cluding recent developments of QCD at non-zero baryon number density. 
We continue with the results obtained so far for the transition tempera- 
ture as well as the temperature dependence of energy and pressure and 
comment on screening and the heavy quark free energies. A major section 
is devoted to the discussion of thermal modifications of hadron properties, 
taking special account of recent progress through the use of the maximum 
entropy method. 
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1 Introduction 

Understanding the properties of elementary particles at high temperature and 
density is one of the major goals of contemporary physics. Through the study 
of properties of elementary particle matter exposed to such extreme conditions 
we hope to learn about the equation of state that controlled the evolution of the 
early universe as well as the structure of compact stars. A large experimental 
program is devoted to the study of hot and dense matter created in ultrarel- 
ativistic heavy ion collisions. Lattice studies of QCD thermodynamics have 
established a theoretical basis for these experiments by providing quantitative 
information on the QCD phase transition, the equation of state and many other 
aspects of QCD thermodynamics. 

Already 20 years ago lattice calculations first demonstrated that a phase 
transition in purely gluonic matter exists ^ |2j and that the equation of state of 
gluonic matter rapidly approaches ideal gas behavior at high temperature 
These observables have been of central interest in numerical studies of the ther- 
modynamics of strongly interacting matter ever since. The formalism explored 
in these studies, its further development and refinement has been presented in 
reviews jl] and the steady improvement of numerical results is regularly pre- 
sented at major conferences 0. 

Rather than discussing the broad spectrum of topics approached in lattice 
studies of QCD thermodynamics we will concentrate here on basic parameters, 
which are of direct importance for the discussion of experimental searches for the 
QCD transition to the high temperature and/or density regime, which generally 
is denoted as the Quark Gluon Plasma (QGP). In our discussion of the QCD 
phase diagram, the transition temperature and the equation of state we will 
also emphasize the recent progress made in lattice studies at non-zero baryon 
number density. A major part of this review, however, is devoted to a discussion 
of thermal modifications of hadron properties, a topic which is of central impor- 
tance for the discussion of experimental signatures that can provide evidence 
for the thermal properties of the QGP as well as those of a dense hadronic gas. 



1.1 QCD Thermodynamics 

A suitable starting point for a discussion of the equilibrium thermodynamics 
of elementary particles interacting only through the strong force is the QCD 
partition function represented in terms of a Euclidean path integral. The grand 
canonical partition function, Z(V, T, fj,f), is given as an integral over the funda- 
mental quark (?/>, ip) an d gluon (A v ) fields. In addition to its dependence on 
volume (V), temperature (T) and a set of chemical potentials (/i/), the parti- 
tion function also depends on the coupling g and on the quark masses to/ for 
/ = 1, .., different quark flavors, 




(1) 
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Here the bosonic fields A v and the Grassmann valued fermion fields ip> "0 obey 
periodic and anti-periodic boundary conditions in Euclidean time, respectively. 
The Euclidean action Se = Sq + Sf contains a purely bosonic contribution (Sq) 
expressed in terms of the field strength tensor, F^ = d ll A u — d u A^ — ig[A^, A v \ , 
and a fermionic part (Sf), which couples the gluon and quark fields through 
the standard minimal substitution, 

Se(V,T, /if) ee S G (V,T) + S F (V,T,» f ) , (2) 

l/T 

S G (V,T) = J dx 4 J d 3 x ^Tr F^F^ , (3) 
o v 

S F (V,T,ii f ) = Jdx 4 jd 3 x ^^(^[^-i^l+wo+m/)^. (4) 

V f=1 

Basic thermodynamic quantities like the pressure (p) and the energy density 
(e) can then easily be obtained from the logarithm of the partition function, 

h = v^3 lnZ ( T '^/) ' (5) 

£ - 3P r _d_ (V_\ , fi x 

T 4 ~ AT VT 4 J|fi X cd M /T ' { > 

Moreover, the phase structure of QCD can be studied by analyzing observ- 
ables which at least in certain limits are suitable order parameters for chiral 
symmetry restoration (rrif — > 0) or deconfinement (my — > oo), i.e. the chiral 
condensate and its derivative the chiral susceptibility, 

&f^ = V^ lnZ{T > V ^ f) ' ^ = |EilnZ(T,V,M/), (7) 
f /=i / 

as well as the expectation value of the trace of the Polyakov loop 1 , 

X 

where the trace is normalized such that Trl = 1. Here L{x) denotes a closed 
line integral over gluon fields which represents a static quark source, 

We may couple these static sources to a constant external field, h, and consider 
its contribution to the QCD partition function. The corresponding susceptibility 
is then given by the second derivative with respect to h, 

XL = V ((L 2 ) - (L) 2 ) , (10) 

X A more formal definition of (L) which leads to a well defined Polyakov loop expectation 
value also in the continuum limit is given in Section 5. 
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where h has been set to zero again after taking the derivatives. 

1.2 Lattice formulation of QCD Thermodynamics 

The path integral appearing in Eq. ^ is regularized by introducing a four di- 
mensional space-time lattice of size N% x N T with a lattice spacing a. Volume 
and temperature are then related to the number of points in space and time 
directions, respectively, 

V = (N a a) 3 , T~ l =N T a , (11) 

and also chemical potentials and quark masses are expressed in units of the 
lattice spacing, jj,f — /i/a, rhf — rrifa. The lattice spacing then does not 
appear explicitly as a parameter in the discretized version of the QCD partition 
function. It is controlled through the bare couplings of the QCD Lagrangian, 
i.e. the gauge coupling 2 g 2 and quark masses rhf. 

At least on the naive level the discretization of the fermion sector is straight- 
forwardly achieved by replacing derivatives by finite differences and by intro- 
ducing dimensionless, Grassmann valued fermion fields. Enforced by the re- 
quirement of gauge invariance the discretization of the gauge sector, however, 
is a bit more involved. Here one introduces link variables U^{x) which are as- 
sociated with the link between two neighboring sites of the lattice and describe 
the parallel transport of the field A from site x to the neighboring site in the /x 
direction x + fia, 

ig J dx^A^x)) , (12) 

where P denotes the path ordering. The link variables U^{x) are elements of 
the SU(3) color group. 

We will not elaborate here any further on details of the lattice formulation 
which is described in excellent textbooks Cj • In recent years much progress 
has also been made in constructing improved discretization schemes for the glu- 
onic as well as fermionic sector of the QCD Lagrangian which greatly reduced 
the systematic errors introduced by the finite lattice cut-off. This improvement 
program and also the improvement of numerical algorithms is reviewed regu- 
larly at lattice conferences and is discussed in review articles [E]. It has been 
crucial also for the calculation of thermodynamic observables and their extrap- 
olation to the continuum limit. The impressive accuracy that can be achieved 
through a systematic analysis of finite cut-off effects on the one hand and the 
use of improved actions on the other hand is apparent in the heavy quark mass 
limit of QCD, i.e. in the SU(3) gauge theory. In particular, for bulk thermo- 
dynamic quantities like the pressure or energy density the discretization errors 
can become huge as these quantities depend on the fourth power of the lattice 



2 In the lattice community it is customary to introduce instead the coupling /3 = 6/g 2 . 
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Figure 1: Discretization errors in the calculation of the pressure of a non- 
interacting gluon gas on lattices with temporal extent N T . Shown are results 
for the standard one plaquette action introduced by K. G. Wilson |U] and the 
renormalization group improved action constructed by Y. Iwasaki ^jj. Both ac- 
tions lead to discretization errors of 0(a 2 ) = D(l/N%). Also shown are results 
obtained with an Symanzik-improved action [111 I12| which has discretization 
errors of C(a 4 ) only. 

spacing a. Nonetheless, improved discretization schemes lead to a large reduc- 
tion of discretization errors and allow a safe extrapolation from lattices with 
small temporal extent N T to the continuum limit (N T — > oo, a — ► 0) at fixed 
temperature T — l/N T a. This is illustrated in Fig.^which shows results for the 
pressure, p(N T ), in a non-interacting gluon gas calculated on a lattice with finite 
temporal extent in comparison to the Stefan-Boltzmann value, psb , obtained in 
the continuum. A similarly strong reduction of cut-off effects can be achieved 
in the fermion sector when using improved fermion actions 0]|B]. 



There are many indications that strongly interacting matter at high temper- 
atures/densities behaves fundamentally different from that at low tempera- 
tures/densities. On the one hand it is expected that the copious production 
of resonances in hadronic matter which will occur in a hot interacting hadron 
gas does set a natural limit to hadronic physics described in terms of ordinary 
hadronic states. On the other hand it is the property of asymptotic freedom 
which suggests that the basic constituents of QCD, quarks and gluons, should 
propagate almost freely at high temperature/densities. This suggests that the 
non-perturbative features characterizing hadronic physics at low energies, con- 
finement and chiral symmetry breaking, get lost when strongly interacting mat- 



2 The QCD phase diagram 
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ter is heated up or compressed. 

Although the early discussions of the phase structure of hadronic matter, 
e.g. based on model equations of state, seemed to suggest that the occurrence 
of a phase transition is a generic feature of strong interaction physics, we know 
now that this is not at all the case. Whether the transition from the low tem- 
perature/density hadronic regime to the high temperature/density regime is 
related to a true singular behavior of the partition function leading to a first or 
second order phase transition or whether it is just a more or less rapid crossover 
phenomenon crucially depends on the parameters of the QCD Lagrangian, i.e. 
the number of light or even massless quark flavors. In particular, whether the 
transition in QCD with values of quark masses as they are realized in nature is 
a true phase transition or not is a detailed quantitative question. The answer 
to it most likely is also dependent on the physical boundary conditions, i.e. 
whether the transition takes place at vanishing or non-vanishing values of net 
baryon number density (chemical potential). 

At vanishing or small values of the chemical potential the crucial control 
parameter is the strange quark mass. Quite general symmetry considerations 
and universality arguments for thermal phase transitions |13| suggest that in 
the limit of light up and down quark masses the transition is first order if also 
the strange quark mass is small enough whereas it is just a continuous crossover 
for strange quark masses larger than a certain critical value, m crlt . This critical 
value also depends on the value of the light quark masses m u , mj. At m crlt the 
transition would be a second order phase transition with well defined universal 
properties, which are those of the three dimensional Ising model |14| . 



3-flavour phase diagram 
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Figure 2: The QCD phase diagram of 3-flavor QCD with degenerate (u,d)-quark 
masses and a strange quark mass m s . 
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The dependence of the QCD phase transition on the number of flavors and 
the values of the quark masses has been analyzed in quite some detail for van- 
ishing values of the chemical potential. The current understanding of this phase 
diagram is summarized in Fig. [3 Its basic features have been established in 
numerical calculations and are in agreement with the general considerations 
based on universality and the symmetries of the QCD Lagrangian. Nonethe- 
less, the numerical values for critical temperatures and critical masses given in 
this figure should just be taken as indicative; not all of them have been deter- 
mined with sufficient accuracy. It is, however, obvious that there is a broad 
range of quark mass values or equivalently pseudo-scalar meson masses, where 
the transition to the high temperature regime is not a phase transition but a 
continuous crossover. The regions of first order transitions at large and small 
quark masses are separated from this crossover regime through lines of second 
order transitions which belong to the universality class of the 3d Ising model 
|15| . The chiral critical line in the small quark mass region has been analyzed 
in quite some detail ^3^3 El- I* 1 the case of three degenerate quark masses 
(3-flavor QCD) it has been verified that the critical point belongs to the Ising 
universality class [TJ]. The critical quark mass m crlt at this point, however, 
is not a universal quantity and is not yet known to a satisfying precision: it 
corresponds to a pseudo-scalar mass varying between about 290 MeV in a stan- 
dard discretization and about 200 MeV in an improved one. This result can 
be extended to the case of non-degenerate quark masses, i.e. m s ^ m; where 
mi = m u = rrid denotes the value of the two light quark masses, which are taken 
to be degenerate. The slope of the chiral critical line in the vicinity of the three 
flavor point can be obtained from a Taylor expansion, 

m c nt _ m CTit + 2 ^Crit _ mj ) . (13) 

This relation has been verified explicitly in a numerical calculation and is consis- 
tent with all existing studies of the chiral transition line |15M16lH7| . A collection 
of results taken from Ref. JHj is shown in Fig. 

If we assume that the linear extrapolation of the chiral critical line, Eq. 1131 
holds down to light quark mass values, which correspond to the physical pion 
mass, we estimate that the ratio of strange to light quark masses at this point is 
(mf lt /rnf ys ) fl= o ~ 5 — 10, depending on the action chosen. This clearly is too 
small to put the physical QCD point, which corresponds to m s /m u ~ 20 into 
the first order regime of the phase diagram. At vanishing chemical potential the 
QCD transition thus most likely is a continuous crossover. 

The phase diagram shown in Fig. for vanishing quark chemical potential 3 
(jti) can be extended to non-zero values of \x. The chiral critical line discussed 
above then is part of a critical surface, m c "*(/i). For small values of the chemical 
potential it can be analyzed using a Taylor expansion of the fermion determinant 
in terms of (/x/T). A preliminary analysis performed for three flavor QCD |18j 



3 The baryon chemical potential is given by /ig = 3/x. 
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Figure 3: The chiral critical line in the light and strange quark mass plane. 
All the results shown have been obtained with a standard staggered fermion 
discretization |18| . 

yields, 

(^) =(^) + 0.21(6)^) 2 + O(( M /T) 4 ) , n/=3. (14) 

The positive slope obtained for m c "*(/i) suggests that at the physical value 
of the strange quark mass a first order transition can occur for values of the 
chemical potential larger than a critical value determined from m crlt (fi) = m s 
(see Fig. EI). 

A first direct determination of the chiral critical point in QCD has been 
performed by Fodor and Katz [TQ. As in the case of a Taylor expansion of 
the QCD partition functions in terms of the chemical potential, which we have 
discussed so far, they have also performed numerical calculations at /i = 0. 
However, they then evaluate the exact ratios of fermion determinants calculated 
at fj, = and \x > and use these in the statistical reweighting of gauge field 
configurations 4 to extend the calculation to fi ^ 0. They find for the chiral 
critical point |H] {v/T) cnt ^ 1.4, hT* = 3^ c "* = (725 ± 35) MeV. Although 
this result still has to be established more firmly through calculations with 
lighter up and down quark masses on larger lattices and with improved actions, 
it also suggests that the dependence of the transition temperature on /j, is rather 
weak. 

An alternative approach to numerical calculations at non-zero values of the 
baryon number density (/i > 0) is based on numerical calculations with imagi- 
nary chemical potentials |22l 1231 123] {hi). This allows straightforward numerical 

4 This approach is well known under the name of Ferrenberg-Swendsen reweighting 1211 and 
finds widespread application in statistical physics as well as in lattice QCD calculations. 
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Figure 4: 3d sketch of the QCD phase diagram ^H| : Shown is the critical surface 
of second order phase transition which separates the regime of first order phase 
transitions at large values of the chemical potential and/or small values of the 
light (m Ui d) and strange (m s ) quark masses from the regime of continuous, 
non-singular transitions (crossover) to the QCD plasma phase. 

calculations for fij > 0. The results obtained in this way then have to be an- 
alytically continued to real values. For small values of the chemical potential 
which have been analyzed so far they turn out to be consistent with the results 
obtained with reweighting techniques. 

Finally we want to note that the discussion of the dependence of the QCD 
phase diagram on the baryon chemical potential in general is a multi-parameter 
problem. As pointed out in Eq. ^ one generally has to deal with independent 
chemical potentials /if for each of the different quark flavors, which control the 
corresponding quark number densities, 



The chemical potentials thus are constrained by boundary conditions which are 
enforced on the quark number densities by a given physical system. In the 
case of dense matter created in a heavy ion collision this is fj, s — due to the 
requirement that the overall strangeness content of the system vanishes. In 
a dense star, on the other hand, weak decays will lead to an equilibration of 
strangeness and it is the charge neutrality of a star which controls the relative 
magnitude of strange and light quark chemical potentials |25| . 

So far we have discussed a particular corner of the QCD phase diagram, 
i.e. the regime of small values of the chemical potential. In fact, the numerical 
techniques used today to simulate QCD with non-vanishing chemical poten- 
tial seem to be reliable for u./T<l and high temperature. Fortunately this is 
the regime, which also seems to be accessible experimentally. On the other 
hand there is the entire regime of low temperature and high density which is 
of great importance in the astrophysical context. In this regime it is expected 
that interesting new phases of dense matter exist. At low temperature and 




z f = 



(15) 
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asymptotically large baryon number densities asymptotic freedom insures that 
the force between quarks will be dominated by one-gluon exchange. As this 
leads to an attractive force among quarks, it seems unavoidable that the naive 
perturbative ground state is unstable at large baryon number density and that 
the formation of a quark-quark condensate leads to a new color-superconducting 
phase of cold dense matter. As this part of the phase diagram is not in the focus 
of our following discussion we refer the interested reader to the many excellent 
reviews which appeared in recent years [23 ■ In Fig. El we just show a sketch of 
the phase diagram of QCD for a realistic quark mass spectrum. The transition 
to the plasma phase is expected to be a continuous crossover (dashed line) for 
small values of the quark chemical potential and turns into a first order tran- 
sition only beyond a critical value for the chemical potential or baryon number 
density. The details of this phase diagram in particular in the low temperature 
regime are, however, largely unexplored in lattice calculations. 



Figure 5: The QCD phase diagram in the T — /i as it might look like in the case 
of non-zero but light up and down quarks and a heavier strange quark. 



3 The transition temperature 

Although a detailed analysis of the QCD phase diagram clearly is of importance 
on its own, it is the thermodynamics in the region of small or vanishing chem- 
ical potential which is of most importance for a discussion of thermodynamic 
conditions created in heavy ion collisions at RHIC or LHC. Current estimates 
of baryon number densities obtained at central rapidities in heavy ion collisions 
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at RHIC PZ] suggest that the baryon chemical potential is below 50 MeV, cor- 
responding to a quark chemical potential of about 15 MeV. We will see in the 
following that at these small values of \i the critical temperature is expected to 
change by less than a percent from that at /i = 0. For this reason, and of course 
also because much more quantitative results are known in this case, we will in 
the following focus our discussion on the thermodynamics at ^ = 0. 

As discussed in the previous section the transition to the high tempera- 
ture phase is continuous and non-singular for a large range of quark masses. 
Nonetheless, for all quark masses the transition proceeds rather rapidly in a 
small temperature interval. A definite transition point thus can be identified, 
for instance through the location of maxima of the susceptibilities of the chiral 
condensate, Xm (Eq. 0), or the Polyakov loop, \l fEci. ITU. While the maxi- 
mum of Xm determines the point of maximal slope in the chiral condensate, \l 
characterizes the change in the long distance behavior of the heavy quark free 
energy (see section EJ). 

On a lattice with temporal extent N T and for a given value of the quark 
mass the susceptibilities define pseudo-critical couplings j3 pc (m q , N T ) which are 
found to coincide within statistical errors. In order to determine the transition 
temperature T c = l/N T a((3 pc ) one has to fix the lattice spacing through the 
calculation of an experimentally or phenomenologically known observable. For 
instance, this can be achieved through the calculation of a hadron mass, m#a, 
or the string tension, era 2 , at zero temperature and the same value of the lattice 
cut-off, i.e. at f3 pc . This yields T c /^/a = l/N T V a a 2 and similarly for a hadron 
mass. In the pure gauge theory the transition temperature has been analyzed 
in great detail and the influence of cut-off effects has been examined through 
calculations on different size lattices and with different actions. From this one 
finds for the critical temperature of the first order phase transition 5 

SU(3) gauge theory : T c /y/a = 0.632 ± 0.002 

T c = (269 ± 1) MeV (16) 

Already the early calculations for the transition temperature in QCD with 
dynamical quark degrees of freedom |31II32| indicated that the inclusion of light 
quarks leads to a significant decrease of the transition temperature. A com- 
pilation of newer results [321 1331 1341 1351 156*| which have been obtained using 
improved lattice regularizations for staggered as well as Wilson fermions is pre- 
sented in Fig. EJ The figure only shows results for 2-flavor QCD obtained from 
calculations with several bare quark mass values. In order to compare calcu- 
lations performed with different actions the results are presented in terms of 
a physical observable, the ratio of pseudo-scalar (pion) and vector (rho) meson 
masses, mps/mv- 

From Fig. it is evident that T c /mv drops with increasing ratio mps/mv, 
i.e. with increasing quark mass. This is not surprising as my, of course, does 

5 This number is a weighted average of the data given in Rcf. 28 29 , including a rephrase- 
ment of r$T c given in Ref. I .' i ( ) I using the string model value a = iz/12 for the 1/R term in the 
heavy quark potential. We also use y/a = 425 MeV for the string tension to set the scale for 
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Figure 6: Transition temperature in 2-flavor QCD in units of my calculated 
on lattices with temporal extent N T — 4 with different fermion actions. The 
large dot drawn for mpg/my — indicates the result of chiral extrapolations 
based on calculations with improved Wilson as well as improved staggered 
|33| fermions. The vertical line shows the location of the physical limit, mps = 
m, = 140 MeV. 



not take on the physical p-meson mass value as long as the quark masses used 
in the calculations are larger than than those realized in nature and the ratio 
m-ps/m-v thus does not attain its physical value (vertical line in Fig. EJ- In 
fact, for large quark masses my will continue to increase with m q while T c will 
remain finite and eventually approach the value calculated in the pure SU (3) 
gauge theory. The ratio T c /my has to approach zero for mps/my = 1. Fig. El 
alone thus does not allow us to quantify the quark mass dependence of T c . 

A simple percolation picture for the QCD transition would suggest that 
T c (m q ) or better T c (mps) will increase with increasing m q ; with increasing m q 
all hadrons will become heavier and it will become more and more difficult to 
excite these heavy hadronic states. It thus becomes more difficult to create a 
sufficiently high particle/energy density in the hadronic phase that can trigger 
a phase (percolation) transition. Such a picture also follows from chiral model 
calculations [23 EH! • ^ n order to check to what extent this physically well mo- 
tivated picture finds support in the actual numerical results obtained in lattice 
calculations of the quark mass dependence of the transition temperature we 
should express T c in units of an observable, which itself is not (or only weakly) 
dependent on m q ; the string tension (or also a hadron mass in the valence quark 
chiral limit 6 ) seems to be suitable for this purpose. In fact, this is what we have 
tacitly assumed when converting the critical temperature of the SU(3) gauge 
theory, T c /^/a ~ 0.63, into physical units as it has been done in Eq. 1161 This 
assumption is supported by the observation that already in the heavy quark 



This often is called the partially quenched limit. 
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Figure 7: The transition temperature in 2 (filled squares) and 3 (circles) flavor 
QCD versus mps/y/o 1 using an improved staggered fermion action (p4- action). 
Also shown are results for 2-flavor QCD obtained with the standard staggered 
fermion action (open squares). The dashed band indicates the uncertainty on 
T c j^fu in the quenched limit. The straight line is the fit given in Eg. 1181 



mass limit the string tension calculated in units of quenched hadron masses, 
e.g. m p /y/a = 1.81 (4) is in good agreement with values required in QCD 
phenomenology, ^/a ~ 425 MeV. 

To quantify the quark mass dependence of the transition temperature one 
may express T c in units of ^fo . This ratio is shown in Fig. as a function 
of mps/^Jcr. As can be seen the transition temperature starts deviating from 
the quenched value for mps < (6 — — 2.5 GeV. We also note that the 

dependence of T c on mpsj \fo is almost linear in the entire mass interval. Such 
a behavior might, in fact, be expected for light quarks in the vicinity of a 2 nd 
order chiral transition where the dependence of the pseudo-critical temperature 
on the mass of the Goldstone-particle follows from the scaling relation 

T c (m») - T c (0) ~ mW* . (17) 

For 2-flavor QCD the critical indices j3 and 5 are expected to belong to the 
universality class of 3-d, 0(4) symmetric spin models and one thus indeed would 
expect to find 2/05 ~ 1.1. However, this clearly cannot be the origin for the 
quasi linear behavior which is observed already for rather large hadron masses 
and, moreover, seems to be independent of n/. In fact, unlike in chiral models 
|23 Ell the dependence of T c on mps turns out to be rather weak. The line 
shown in Fig. is a fit to the 3- flavor data, 

(%) =(%) +0 -^( !! w) • ( 18 ) 



It seems that the transition temperature does not react strongly to changes 
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of the lightest hadron masses. This favors the interpretation that the contribu- 
tions of heavy resonance masses are equally important for the occurrence of the 
transition. In fact, this also can explain why the transition still sets in at quite 
low temperatures even though all hadron masses, including the pseudo-scalars, 
attain masses of the order of 1 GeV or more. Such an interpretation also is 
consistent with the weak quark mass dependence of the critical energy density 
which one finds from the analysis of the QCD equation of state as we will discuss 
in the next section. 

In Fig. we have included results from calculations with 2 and 3 degenerate 
quark flavors. So far such calculations have mainly been performed with stag- 
gered fermions. In this case also a simulation with non-degenerate quarks (a 
pair of light u,d quarks and a heavier strange quark) has been performed. Un- 
fortunately, the light quarks in this calculation are still too heavy to represent 
the physical ratio of light u,d quark masses and a heavier strange quark mass. 
Nonetheless, the results obtained so far suggest that the transition tempera- 
ture in (2+l)-flavor QCD is close to that of 2-flavor QCD. The 3-flavor theory, 
on the other hand, leads to consistently smaller values of the critical temper- 
ature, T c (n,f = 2) — T c (rif — 3) ~ 20 MeV. Extrapolations of the transition 
temperatures to the chiral limit gave 

clover-improved Wilson 
fermions |3*3~] 
improved staggered 
fermions |35| 
improved staggered 
fermions |35j 

Here m p has been used to set the scale for T c . Although the agreement between 
results obtained with Wilson and staggered fermions is striking, one should 
bear in mind that all these results have been obtained on lattices with temporal 
extent N T = 4, i.e. at rather large lattice spacing, a ~ 0.3 fm. Moreover, there 
are uncertainties involved in the ansatz used to extrapolate to the chiral limit. 
We thus estimate that the systematic error on the value of T c /m p still is of 
similar magnitude as the purely statistical error quoted above. 

As mentioned already in the previous section, first studies of the dependence 
of the transition temperature on the chemical potential have been performed re- 
cently using either a statistical reweighting technique ^1121310131 to extrapolate 
from numerical simulations performed at = to > or performing simula- 
tions with an imaginary chemical potential [231 124| /x/ the results of which are 
then analytically continued to real fi. To leading order in a 2 one finds 

Tc(v) = ( 1 - 0.0056(4)(^ B /T) 2 Ref. (imaginary n) ( 
T c (0) ~ 1 1 - 0.0078(38) {^b/T? Ref. 40^ (C(^ 2 ) reweighting) ' [ ' 

These results are consistent with the (2+1) -flavor calculation performed with an 
exact reweighting algorithm |19M20| . The result obtained for T c (fi) in this latter 
approach is shown in Fig. |S1 The dependence of T c on the chemical potential 



2 - flavor QCD : T c = 



(171 ± 4) MeV, 
(173 ± 8) MeV, 



3 - flavor QCD : T c = (154 ± 8) MeV, 
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Figure 8: The /i-dependence of the transition temperature for (2+l)-flavor QCD 



is rather weak. We stress, however, that these calculations have not yet been 
performed with sufficiently light up and down quark masses and a detailed 
analysis of the quark mass dependence has not yet been performed. The [i- 
dependence of T c (/z)/T c (0) is expected to become stronger with decreasing quark 
masses (and, of course, vanishes in the limit of infinite quark masses). 

4 The equation of state 

One of the central goals in studies of the thermodynamics of QCD is, of course, 
the calculation of basic thermodynamic quantities and their temperature de- 
pendence. In particular, one wants to know the pressure and energy density 
which are of fundamental importance when discussing experimental studies of 
dense matter. Besides, they allow a detailed comparison of different computa- 
tional schemes, e.g. numerical lattice calculations and analytic approaches in 
the continuum. 

At high temperature one generally expects that due to asymptotic free- 
dom these observables show ideal gas behavior and thus are directly propor- 
tional to the basic degrees of freedom contributing to thermal properties of 
the plasma, e.g. the asymptotic behavior of the pressure will be given by the 
Stcfan-Boltzmann law 



Perturbative calculations 01] of corrections to this asymptotic behavior are, 
however, badly convergent and suggest that a purely perturbative treatment 
of bulk thermodynamics is trustworthy only at extremely high temperatures, 



123- 



^ = ( 16+10 -^4) 



(20) 
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i.e. several orders of magnitude larger than the transition temperature to the 
plasma phase. In analytic approaches one thus has to go beyond perturbation 
theory which currently is being attempted by either using hard thermal loop 
resummation techniques in combination with a variational ansatz |42l I4H| or 
perturbative dimensional reduction combined with numerical simulations of the 
resulting effective 3-dimensional theory |441 14"5] . 

The lattice calculation of pressure and energy density is based on the stan- 
dard thermodynamic relations given in Eq. For vanishing chemical potential 
the free energy density is directly given by the pressure, / = —p. As the par- 
tition function itself is not directly accessible in a Monte Carlo calculation one 
first takes a suitable derivative of the partition function, which yields a calcula- 
ble expectation value, e.g. the gauge action. After renormalizing this observable 
by subtracting the zero temperature contribution it can be integrated again to 
obtain the difference of free energy densities at two temperatures, 



The lower integration limit T is chosen at low temperatures so that p/T 4 is small 
and may be ignored. This easily can be achieved in an SU(3) gauge theory where 
the only relevant degrees of freedom at low temperature are glueballs. Even 
the lightest ones calculated on the lattice have large masses, m& ~ 1.5 GeV. 
The free energy density thus is exponentially suppressed up to temperatures 
close to T c . In QCD with light quarks, however, the dominant contribution 
to the free energy density comes from pions. In the small quark mass limit 
also T Q has to be shifted to rather small temperatures. At present, however, 
numerical calculations are performed with rather heavy quarks and also the 
pion contribution thus is strongly suppressed below T c . 

In Fig. we show results for the pressure obtained in calculations with 
different numbers of flavors |46| . At high temperature the magnitude oip/T 
clearly reflects the change in the number of light degrees of freedom present in 
the ideal gas limit. When we rescale the pressure by the corresponding ideal gas 
values it becomes, however, apparent that the overall pattern of the temperature 
dependence of p/T 4 is quite similar in all cases. The figure also shows that the 
transition region shifts to smaller temperatures as the number of degrees of 
freedom is increased. As pointed out in the previous section such a conclusion, 
of course, requires the determination of a temperature scale that is common to 
all QCD -like theories which have a particle content different from that realized 
in nature. We have determined this temperature scale by assuming that the 
string tension is flavor and quark mass independent. 

Other thermodynamic observables can be obtained from the pressure using 
suitable derivatives. In particular one finds for the energy density, 




(21) 




(22) 
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Figure 9: The pressure in QCD with two and three degenerate quark flavors as 
well as two light and a heavier (strange) quark ^jg\. For nj ^ calculations 
have been performed on an N T = 4 lattice using improved gauge and staggered 
fermion actions with a quark mass m/T = 0.4. Cut-off effects in these calcu- 
lations are expected to be of the order of 20% and the pressure is expected to 
become correspondingly larger once a proper extrapolation to the continuum 
limit can be performed on larger lattices. At high temperature the influence of 
a non-zero quark mass is expected to be small. In the case of the SU(3) pure 
gauge theory the continuum extrapolated result is shown. 

In Fig. ^| we show results for the energy density 7 obtained from calculations 
with staggered fermions and different number of flavors. Unlike the pressure the 
energy density rises rapidly at the transition temperature. Although the results 
shown in this figure correspond to quark mass values in the crossover region 
of the QCD phase diagram the transition clearly proceeds rather rapidly. This 
has, for instance, also consequences for the velocity of sound, = dp/de, which 
becomes rather small close to T c . The velocity of sound is shown in Fig. El The 
comparison of results obtained from calculations in the SU (3) gauge theory [2Hj 
and results obtained in simulations of 2-flavor QCD using Wilson fermions |4*7j 
with different values of the quark mass shows that the temperature dependence 
of v s is almost independent of the value of the quark mass. 

Also shown in Fig. ^| is an estimate of the critical energy density at which 
the transition to the plasma phase sets in. In units of T% the transition takes 
place at e/Tjf ~ 6 ± 2 which should be compared with the corresponding value 
in the SU(3) gauge theory [2Hj, ^/T^ ~ 0.5. Although these numbers differ 
by an order of magnitude it is rather remarkable that the transition densities 



7 This figure is based on data from Ref. 1461 obtained for bare quark masses m/T = 0.4. 
The energy density shown does not contain a contribution which is proportional to the quark 
mass and thus vanishes in the chiral limit. 
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Figure 10: The energy density in QCD with 2 and 3 degenerate quark flavors. 
Also shown is a sketch of the expected form of the energy density for QCD with 
a fixed strange quark mass m s ~ T c (see also remarks on cut-off effects in the 
caption of Fig. The arrows indicating the energy densities reached in the 
initial stage of heavy ion collisions at the SPS, RHIC and in the future also at 
the LHC are based on the Bjorken formula |48j . 



expressed in physical units are quite similar in both cases; when moving from 
large to small quark masses the increase in e/T 4 is compensated by the decrease 
in T c . This result thus suggests that the transition to the QGP is controlled by 
the energy density, i.e. the transition seems to occur when the thermal system 
reaches a certain "critical" energy density. In fact, this assumption has been 
used in the past to construct the phase boundary of the QCD phase transition 
in the T — [i plane. 

Also at non-vanishing baryon number density, the pressure as well as the 
energy density can be calculated along the same line outlined above by us- 
ing the basic thermodynamic relation given in Eq. [H] Although the statisti- 
cal errors are still large, a first calculation of the /x-dependence of the transi- 
tion line indeed suggests that e(T c (/i),/i) varies only little with increasing /i, 
e(T c (/x),/i) - e(T c (0),0) = (1.0 ± 2.2)^T C 2 (0) 40;. First calculations of the in- 
dependence of the pressure in a wider temperature range have recently been 
performed using the reweighting approach for the standard staggered fermion 
formulation |49j as well as the Taylor expansion for an improved staggered 
fermion action up to 0((/z/T) 4 ) jSUj. This shows that the behavior of bulk 
thermodynamic observables follow a similar pattern as in the case of vanishing 
chemical potential. For instance, the additional contribution to the pressure, 
Ap/T 4 = (p/T 4 )^/ T — (p/T 4 ) Al=0 rapidly rises at T c and shows only little tem- 
perature variation for T/T c >l.b. In this temperature regime the dominant con- 
tribution to the pressure arises from the contribution proportional to ([i/T) 2 
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Figure 11: The velocity of sound in the SU(3) gauge theory [2Hj and in 2-fLavor 
QCD 47 . In the latter case we show results from calculations with Wilson 
fcrmions performed at different values of the quark mass. 



which also is the dominant contribution in the ideal gas limit as long as fi/T<l, 

(£L-(£L=T(f) 2 +^)* ■ <»> 

Here (p/T 4 )^^ is given by Eq. [2UJ It also turns out that at non-zero values 
of the chemical potential the cut-off effects in bulk thermodynamic observables 
are of the same size as at /i = 0. Further detailed studies of the behavior of the 
pressure and the energy density thus will require a careful extrapolation to the 
continuum limit and/or the use of improved gauge and fermion actions. 

5 Heavy quark free energies 

Heavy quark free energies play a central role in our understanding of the QCD 
phase transition, the properties of the plasma phase and the temperature depen- 
dence of the heavy quark potential. The heavy quark free energies £Q are defined 
as the QCD partition functions of thermal systems containing n (n) static quark 
(anti-quark) sources located at positions x = {xi}" =1 and x = {xi}™ =1 , 

Zfrfi)(y,T,x,x) = exp(-F( ,l <")(y,T,x,x)/T) (24) 

/n n 
VA v V^Vi\> e SB(v.T. Mf ) -Q Tr L ^ -q Tr L t(^.) ) 

where the Polyakov Loop, L(x), has been defined in Eq. |pj The expectation 
value of the product of Polyakov loops gives the difference in free energy due to 
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the presence of static gg-sources in a thermal heat bath of quarks and gluons, 

nTrL(f 4 )n TrLt (^)/ = z\ ViT) ( 25 ) 

= cxp(-A J F( n ' fi )(V r ,r,x,x)/T) 

ee exp(-(^ R )(^T,x,x)-F(F,T))/T) , 

where Z(V,T) is the QCD partition function defined in Eq. ^ In particular, 
one considers the two point correlation function (n = n = 1), 

G i (|f-yl,T) = (TrL(f)TrLt(y)) , (26) 

and the Polyakov loop expectation value, which can be defined through the large 
distance behavior of Gl, 



(L) = lim y/G L (r,T) , r = \x-y\ . (27) 

r — >oo 

These observables elucidate the deconfining features of the transition to the high 
temperature phase of QCD. 



5.1 Deconfinement order parameter 

The Polyakov loop expectation value has been introduced in Section 1 as an 
order parameter for deconfinement in the heavy quark mass limit of QCD, i.e. 
the SU(3) gauge theory. Like in statistical models, e.g. the Ising model, it is 
sensitive to the spontaneous breaking of a global symmetry of the theory under 
consideration. In the case of the SU(3) gauge theory this is the global Z(3) 
centre symmetry |51j . In the presence of light dynamical quarks this symmetry 
is explicitly broken and in a strict sense the Polyakov loop looses its property as 
an order parameter. Through its relation to the two point correlation function, 
Eq. E| it however still is the free energy of a static quark placed in a thermal 
heat bath, 

F q (T) = -Tln«L» . (28) 

In the low temperature, confined regime (L) is small and the free energy thus 
is large. It is infinite only for an SU(3) gauge theory, i.e. for QCD in the limit 
m q — > oo. In the high temperature regime however, (L) becomes large and F q 
decreases rapidly when crossing the transition region. 

The static quark sources introduced in the QCD partition function through 
the line integral defined in Eq. also introduce additional ultraviolet diver- 
gences, which require a proper renormalization. For the lattice regularized 
Polyakov loop this can be achieved through a renormalization of the tempo- 
ral gauge link variables Uofaojri), i.e. 



L H = ]JZ L (g 2 )U (i,n) 

i=l 



(29) 
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The renormalization constant Zr 1 (g 2 ) can, for instance, be determined by 
normalizing the two point correlation functions such that the resulting free en- 
ergy of a color singlet quark anti-quark pair at short distances coincides with 
the zero temperature heavy quark potential. This also insures that divergent 
self energy contributions to the Polyakov loop expectation value, defined by 
Eq. 1271 get removed and that the heavy quark free energy can unambiguously 
be defined also in the continuum limit |52| . 

We stress that it is conceptually appealing to define the renormalization con- 
stant for the Polyakov loop in terms of color singlet free energies. Nonetheless, 
this can also be achieved through the gauge invariant two point correlation func- 
tions, Gl(t,T), which define so called color averaged free energies. A qq-pair 
placed in a thermal heat bath cannot maintain its relative color orientation. The 
entire thermal system (gg-pair + heat bath) will be colorless and the gg-pair 
can change its orientation in color space when interacting with gluons of the 
thermal bath. The Polyakov loop correlation function thus has to be considered 
as a superposition of contributions arising from color singlet (F\ ) and color octet 
(Fg) contributions to the free energy £Q, 

e -F(W)(r,T)/T = 1 e -F!(r,r)/T + ^ e -F 8 (r,T)/T _ ^ 

9 9 

At short distances the repulsive octet term is exponentially suppressed and the 
contribution from the attractive singlet channel will dominate the heavy quark 
free energy, 

FW(r,T) = Fx^T) _ 
T T 11 

Q 2 1 

= — - - ln9 for rT«l , (31) 
37r rl 

where the last equality gives the perturbative result obtained from 1-gluon ex- 
change at short distances. When normalizing the color singlet free energy at 
short distances such that it coincides with the zero temperature heavy quark 
potential Vq q {r) the corresponding color averaged free energy thus will differ by 
an additive constant, 

\im{F^' l \r,T) - F x {r,T)) =Tln9 for all T. (32) 

Using this normalization condition the renormalized Polyakov loop order param- 
eter has been determined for the SU(3) gauge theory. It is shown in Fig.^J As 
the deconfmement phase transition is first order in the SU(3) gauge theory the 
order parameter is discontinuous at T c . From the discontinuity, (L)(T C + ) ~ 0.4, 
one finds for the change in free energy F 00 (T C + ) ~ 0.9T C ~ 250 MeV. 

In QCD with light quarks the renormalization program outlined above has 
not yet been performed in such detail as practically all studies of the heavy 
quark free energy have been performed on rather coarse lattices with a small 
temporal extent, N T = 4. Nonetheless, normalizing the free energies obtained 
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Figure 12: Temperature dependence of the renormalized Polyakov loop expec- 
tation value in the pure SU (3) gauge theory determined from the asymptotic 
behavior of color singlet free energies on lattices of size 32 3 x N T . 



in such calculations at the shortest distance presently available (rT — 1/N T ~ 
0.25) to the zero temperature Cornell potential does seem to be a reasonable 
approximation [^| (see Fig. 113(1 . Also in this case the free energy at T c takes 
on a similar value as in the pure gauge theory. 



5.2 Heavy quark potential 

The change in free energy due to the presence of a static quark anti-quark pair 
is given by the two point correlation function defined in Eq. 1261 In the zero 
temperature limit the free energies, F^(r,T), determined from G L (r,T) de- 
fine the heavy quark potential. Also at non zero temperature the free energies 
exhibit properties expected from phenomcnological discussions of thermal mod- 
ifications of the heavy quark potential. In the pure gauge theory m — > oo) 
the free energies diverge linearly at large distances in the low temperature, 
confinement phase. The coefficient of the linear term, the string tension at 
finite temperature, decreases with increasing temperature and vanishes above 
T c . In the deconfined phase the free energies exhibit the behavior of a screened 
potential. At large distances they approach a constant value i 7, 00 (T') at an ex- 
ponential rate. This exponential approach defines a thermal screening mass. 
For finite quark masses the free energies show the expected string breaking be- 
havior; at large distances they approach a finite value at all temperatures, i.e. 
also below T c . This asymptotic value rapidly decreases with decreasing quark 
mass and increasing temperature v 35, as can be seen in Fig. 1131 where we show 
.fVi 1 )^ T) in QCD with three light quark degrees of freedom [SSj. In the lower 
part of Fig. 1131 we show the change in free energy needed to separate a quark 
anti-quark pair from a distance typical for the radius of a cc bound state, i.e. 
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Figure 13: The color averaged heavy quark free energy (V = F at various 
temperatures in the low temperature phase of three flavor QCD (upper figure). 
The free energy has been normalized at the shortest distance available (r = 1/T) 
to the zero temperature Cornell potential, V{r)/^fa = — a/r^/a + r^/a with 
a = 0.25±0.05 (solid band). The lower figure shows the temperature dependence 
of the change in free energy when separating the quark anti-quark pair from a 
distance r = 0.5/ \fa to infinity. To set the scale we use ^/a = 425 MeV. 



r = 0.5/^/cr ~ 0.2 fm, to infinity. 

The change in free energy induced by a heavy quark anti-quark pair often 
is taken to be the heavy quark potential at finite temperature. Of course, this 
is not quite correct and care has to be taken when using the free energies in 
phenomenological discussions of thermal properties of heavy quark bound states. 
In this case we would like to know the energy needed to break up a color singlet 
state formed by a gq-pair. As pointed out in the previous subsection, F^ 1 ' 1 ^ 
represents a color averaged free energy. It may be decomposed in terms of 
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the corresponding singlet and octet contributions, which has indeed been done 
for 5/7(2) and SU(S) gauge theories ESI- However, even then one has 
to bear in mind that the free energy is the difference between an energy and 
entropy contributions, F = E — TS. Of course, from a detailed knowledge of 
the temperature dependence of the free energy at fixed separation of the qq-p&vc 
we can in principle determine the entropy and energy contributions separately, 

dF(r,T) dF{r, T)/T 
S = qjT- , E--T — . (33) 

Such an analysis does, however, not yet exist. 



6 Thermal modifications of hadron properties 

6.1 QCD phase transition and the hadron spectrum 

The non-perturbative structure of the QCD vacuum, in particular confinement 
and chiral symmetry breaking, determine many qualitative aspects of the hadron 
mass spectrum. Also the actual values of light and heavy quark bound state 
masses depend on the values of the chiral condensate and the string tension, 
respectively. As these quantities change with temperature and will change dras- 
tically close to the QCD transition temperature it also is expected that the 
properties of hadrons, e.g. their masses and widths, undergo drastic changes at 
finite temperature. 

Investigations into the nature of hadronic excitations are interesting for var- 
ious reasons in the different temperature regimes. Below T c temperature de- 
pendent modifications of hadron masses and widths may lead to observable 
consequences in heavy ion collision experiments e.g. a pre-deconfined dilep- 
ton enhancement due to broadening of the p-resonance or a shift of its mass 
|5fi| . At temperatures around the transition temperature the (approach to the) 
restoration of chiral symmetry should reflect itself in degeneracies of the hadron 
spectrum. First evidence for this has, indeed, been found early in lattice cal- 
culations of hadronic correlation functions |57|. In the plasma phase the very 
nature of hadronic excitations is a question of interest. Asymptotic freedom 
leads one to expect that the plasma consists of a gas of almost free quarks and 
gluons. The lattice results on the equation of state, however, have shown al- 
ready that this is not yet the case for the interesting temperature region. While 
quasi-particle models and HTL resummed perturbation theory are able to re- 
produce the deviations from the ideal gas behaviour observed in the equation 
of state at temperatures a few times T c it remains to be seen whether they also 
account properly for hadronic excitations. This question arises in particular as 
the generally assumed separation of scales 1/T <C l/(gT) <C l/(g 2 T) does not 
hold for temperatures quite a few times T c . 

In the previous section we have discussed modifications of the heavy quark 
free energy which indicate drastic changes of the heavy quark potential in the 
QCD plasma phase. As a consequence, depending on the quark mass, heavy 
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quark bound states cannot form above certain "critical" temperatures [SB . Sim- 
ilarly it is expected that the QCD plasma cannot support the formation of light 
quark bound states. In the pseudo-scalar sector the disappearance of the light 
pions clearly is related to the vanishing of the chiral condensate at T c . For 
T > T c the pions would be no longer (nearly massless) Goldstone bosons. In 
the plasma phase one thus may expect to find only massive quasi-particle ex- 
citations in the pseudo-scalar quantum number channel. However, also below 
T c it is expected that the gradual disappearance of the spontaneous breaking of 
chiral flavor symmetry as well as the gradual effective restoration of the axial 
Ua (1) symmetry may lead to thermal modifications of hadron properties. While 
the breaking of the SUL{rif) x SUn(nf) flavor symmetry leads, for instance, to 
the splitting of scalar and pseudo-scalar particle masses, the C^t(l) symmetry 
breaking is visible in the splitting of the pion and the rj meson. 

More general, thermal modifications of the hadron spectrum should be dis- 
cussed in terms of modifications of hadronic spectral functions which describe 
the thermal average over transition matrix elements between energy eigenstates 
(E n ) with fixed quantum numbers (H) |59| . 

vh(u,P, T) = ^r J2 e- E "^ /T (l - e-/ T ) 5{u: + E n (p) E rn (p)) ■ 

* ' n,rn 

|(n|J ff (0)|m)| 2 . (34) 

These spectral functions in turn determine the structure of Euclidean correla- 
tion functions, Gh{t,x). Numerical studies of Gh{t,x) at finite temperature 
thus will allow to learn about thermal modifications of the hadron spectrum, al- 
though in practice it is difficult to reconstruct the spectral functions themselves. 
Here recent progress has been reached through the application of the maximum 
entropy method (MEM). Before describing these developments and presenting 
recent lattice results we will start in the next subsection with a presentation 
of some basic field theoretic background on hadronic correlation functions and 
their spectral representation. 

6.2 Spatial and temporal correlation functions, hadronic 
susceptibilities 

6.2.1 Basic field theoretic background 

Numerical calculations of hadronic correlation functions are carried out on Eu- 
clidean lattices i.e. one uses the imaginary time formalism. This holds also for 
zero temperature computations in which case the T — > limit has to be taken. 
The formalism has been worked out in detail in textbooks 021 HSU , however, 
for readers' convenience we have collected some formulae in the Appendix. 
Hadronic correlation functions in coordinate space, Gh(t, x), are defined as 



Gh(t,x) = (J h (t, x)4(0,0)) 



(35) 
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where the hadronic, e.g. mesonic, currents Jh(t,x) = q(r,x)T hq{t,x) contain 
an appropriate combination of 7- matrices, Th, which fixes the quantum num- 
bers of a meson channel; i.e., T H = 1, 75, 7^, 7^75 for scalar, pseudo-scalar, 
vector and pseudo-vector channels, respectively. 

On the lattice, at zero temperature, one usually studies the temporal corre- 
lator at fixed momentum p. Save possible subtractions, the correlator is related 
to the spectral function (see Appendix), (Th(po,P), by means of 

r+00 

g h( t ,P) = / dp a H (p ,p)K(p ,T) , (36) 
Jo 

where the kernel 

cosher - 1/2T)] 
K{P0 ' T)= sinh( P0 /2T) (3?) 

describes the propagation of a single free boson of mass M = po. At zero tem- 
perature, for large temporal separations the correlation function is dominated 
by the exponential decay due to the lightest contribution to the spectral function 
in a given channel. 

At finite temperature, studies of the temporal correlator are hampered by 
the limited extent of the system in the time direction. Therefore most (lattice) 
analyses have been concentrating on the spatial correlation functions. These 
depend of course on the same (temperature dependent) spectral density but are 
different Fourier transforms of it. Projecting onto vanishing transverse momen- 
tum and vanishing Matsubara frequency one obtains 

= r ^ r dpo^^p.) . (38) 



2?T 7-oo Po 

In addition it is quite common in lattice calculations to analyze hadronic sus- 
ceptibilities which are given by the space-time integral over the Euclidean cor- 
relation functions, 

XH = I ' dr Gfj (t, 0) . (39) 
Jo 

The susceptibilities have a particularly simple relation to the spectral function, 

Xh = 2 / dp . (40) 

Jo Po 

Unfortunately, these susceptibilities are generally ultraviolet divergent and the 
above defined integrals should be cut-off at some short distance scale. Rather 
than doing this one can consider a closely related quantity, which provides a 
smooth, exponential cut-off for the ultraviolet part of the spectral function and 
is given by the thermal correlation function at tT = 1/2, 

OUmSl-f ■ (41) 
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In the case of a free stable boson of mass Mjj the spectral function is a pole 
<r H (po,P) = |(0| Jh \H(p))\ 2 e(p )6(p 2 -p 2 - M%) . (42) 

Correspondingly, the imaginary time correlator, projected to vanishing momen- 
tum, decreases with the mass (modulo periodicity) as 

pT,* 1 cosh[M g (r - 1/2T)} 

Likewise, in this case the spatial correlator is also decaying with the mass 

G s H (z)~^-eM~M H z) . (44) 

In this simple case the exponential fall-off of spatial and temporal correlation 
functions thus carry the same information on the particle mass. Interactions in 
a thermal medium, however, are likely to alter the dispersion relation to 

u 2 (p, T) = M 2 H + f + n(p, T) (45) 

with the temperature dependent vacuum polarization H(p,T). In the simplest 
case, one can perhaps assume that the temperature effects can be absorbed 
into a temperature dependent mass Mg-(T) and a coefficient A(T) which might 
also be temperature dependent and different from I, 

u? (p, T) ~ Ml (T) + A 2 (T)p 2 (46) 

In this case, at zero momentum the temporal correlator will decay with the 
so-called pole mass Mh(T), 

G T H (T,6)~eM-M H (T)T) , (47) 

whereas the spatial correlation function has an exponential fall-off, 

G s H (z)~exp(-M s H c (T)z) , (48) 

determined by the screening mass Mf£(T) = Mh(T)/A(T) which differs from 
the pole mass if A(T) ^ 1. In this simple case also the susceptibility defined in 
Eq. as well as the central value of the temporal correlation function defined 
in Eq. ^]are closely related to the pole mass, 

I 

XH 



MKT) ' 

G£(r = l/2T,0) - — — - — - , (49) 



respectively 8 



8 Note that for the spatial correlation functions the relevant symmetry to classify states 
no longer consists of the SO(3) group of rotations but rather is an SO(2) X Z2 due to the 
asymmetry between spatial and temporal directions. This leads to non-degeneracies between 
e.g. po and p x , y excitations 1021 . 
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(a) (b) 

o 



Figure 14: The self-energy diagrams for free quarks (a) and in the HTL approx- 
imation (b). 

Table 1: Coefficients an and bn for the free continuum as well as coefficients 
a}^ and d 1 ^ for the free lattice correlation function, Eqs. 1511581 in the various 
channels, d is given by A/ sinli £. 
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The opposite limit to the case of a free stable boson is reached for two freely 
propagating quarks contributing to the spectral density. Here, to leading order 
perturbation theory the evaluation of the meson correlation function amounts 
to the evaluation of the self-energy diagram shown in Fig. 114b in which the 
internal quark lines represent bare quark propagators Sf(iu n ,p) 03 . 

For massless quarks the spectral density in the mesonic channel H is then 
computed from Eg. I9UI as. 



&h{po,P) = (Po - P 2 ) a H 




(50) 



where p = \f^P and an depends on the channel analyzed (see Tabled for some 
selected values). In the limit of vanishing momentum p the spectral density is 
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also known [54 for quarks with non-vanishing masses to, 



<Jh(po-,0) = £| vl e(Po " 4m 2 ) tanh j i I 




(51) 



The coefficients are also given in Tabled For massless quarks closed analytic 
expressions can be given for both, the temporal as well as the spatial correlator, 
e.g. for the pion 

G T lfr ee (T) 5)/T 3 = ^ _ 2tT) 1±J^IA + 2N c ^Ill ( 52 ) 



sin J (27rTr) 



~sin z (2nTT) 



G 



S,frec 



A r T 



[1 + 2nTz coth(2vrTz)] 



4ttz 2 sinh(27rTz) 
e - M «™" a with AfL ree = 2irT 



(53) 



In this free field limit the susceptibility X-f/> defined in Eo. 1391 is divergent, while 
the temporal correlator at tT — 1/2, of course, stays finite 



G 



T.frce 
H 



(1/2,0) /T 3 = a H N c /3 



(54) 



It, however, is neither related to the screening mass nor does this finite value 
have anything to do with the existence of a pole mass. 

The above discussion can be extended to the leading order hard thermal 
loop (HTL) approximation |64| using dressed quark propagators and vertices 
for the calculation of the self-energy diagrams as indicated in Fig. 114b . The 
corresponding quark spectral functions are given in the Appendix. 

In the interacting case it has been argued |65| from dimcnsionally reduced 
QCD that the simple relation between screening mass and lowest Matsubara 
frequency will be changed, to leading order, into 



M sc (T) = 2nT- 



2 - In 



C F g 2 (T) 



4tt 2 



(55) 



Further corrections arc expected to arise from non-perturbativc effects, e.g. from 
the fact that space-like Wilson loops obey an area law Correspondingly, 
the logarithmic potential used to arrive at the estimate given in Eg. 1551 will turn 
into a linear rising one at large distances. The effect of this is expected to be 
marginal on the screening mass, however, will be large for so-called spatial wave 
functions. 



6.2.2 Lattice results 



As has been mentioned already in the previous section, the temporal separa- 
tion of the operators is limited by the inverse temperature and the traditional 
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Figure 15: Screening masses in quenched QCD calculated with Wilson 

and staggered fermions. Below T c the masses are normalized to T c , above to T. 



approach of extracting a ground state mass from the long distance behavior 
of a Euclidean correlation function thus is not possible. Many lattice studies 
therefore have concentrated on the analysis of spatial correlation functions. As 
outlined in the previous section the exponential decrease in the spatial distance 
then defines a screening mass, which in certain limiting cases indeed may be 
related to the pole mass. 

Fig. 1151 shows one representative set of lattice results on screening masses 
calculated in the quenched approximation and using the Wilson as well as the 
staggered discretization scheme for the fermion action. So far, below T c the 
screening masses themselves have not shown any significant temperature de- 
pendence. This holds for quenched staggered |H7] and Wilson |61M68llr?51l70| as 
well as for dynamical quarks [?7J E2 Thus, there also is no substantial differ- 
ence between the screening and the zero temperature masses. This seems to be 
in accord with sum rule predictions |73j but in conflict with chiral perturbation 
theory [77]. It has been argued [75] that the masses should have the same tem- 
perature dependence as the chiral condensate. However, also the lattice results 
for this quantity when normalized to its zero temperature value do not exhibit 
marked deviations from unity up to temperatures close to T c , see Fig. 1161 Here it 
may be of particular importance to perform numerical calculations with lighter 
dynamical quark masses in order to get better control over the influence of light 
virtual quarks on the chiral condensate and the resonance decay of hadronic 
bound states. 

At temperatures (slightly) above T c the temporal and spatial correlation 
functions of chiral partners reflect the restoration of chiral SU^nj) x SUfi(nf) 
symmetry. In particular, the vector and axial vector channel become degenerate 
[571 1571 16% 1571 1771 17T1 [77117^ independent of the discretization 

and of the number of dynamical flavors being simulated. The same degeneracy is 
observed within errors in the pseudo-scalar and isoscalar scalar (a/fo) channel 
although the latter, for technical reasons, is difficult to access on the lattice. 
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Figure 16: The chiral condensate normalized to its zero temperature value for 
a variety of flavor numbers. 



Nevertheless, screening masses obtained from fits to correlation functions [52*] 
and susceptibilities (at finite lattice spacing) 83, 84 lead to a consistent picture. 

A more intricate question is the one concerning the degeneracy of the pseudo- 
scalar and scalar isovector (5 /do) channel [85|. In three-flavor QCD this degen- 
eracy follows already from the above mentioned chiral one. For two flavors in 
the chiral limit it detects the effective restoration of the anomalous Ua(1). Al- 
though the Ua{\) is explicitly broken by perturbative quantum effects it might 
become effectively restored non-perturbatively if topologically non-trivial zero- 
modes of the Dirac-matrix are absent 9 . If this happens already at temperatures 
below or at T c the chiral transition for two flavors will even be first order |13|. 
Getting control over the zero-modes in lattice calculations is, however, compli- 
cated. They are particularly sensitive to discretization effects and the continuum 
as well as the chiral limit have to be controlled. The use of improved actions 
with better chiral properties at finite lattice spacing thus is important and the 
most convincing evidence so far has therefore been obtained in calculations with 
the domain wall fermion discretization |86j (see Fig. I17fl . This study suggests 
that the £7a(1) symmetry is not yet restored effectively at the transition tem- 
perature, a conclusion which was cautiously drawn also from earlier attempts 
utilizing standard staggered actions [71 1521 1551 151 187] . 

At high temperature, the screening masses are expected to approach the 
free quark propagation limit, Eg. 1531 In fact, already at temperatures as low as 
about 1.5 T c , the results for the vector channel are close to this value, M sc = 
2irT . Nonetheless, Fig. ["""5] still indicates some 10 - 15 % deviations from free 
quark behavior which are of similar size than e.g. in the equation of state. As 

9 In the quenched case, contributions from topologically non-trivial gauge configurations 
are not suppressed by powers of the quark mass arising from the fermion determinant 88 . 
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Figure 17: The difference of the ir and 5 susceptibility in two-flavor QCD as 
function of the quark mass in lattice units. The data |8fi| were obtained at two 
values of the coupling constant which correspond to temperatures slightly above 
T c . The lines are fits of the form cq + c^mj which is the quark mass dependence 
expected in the continuum. A vanishing intercept in the chiral limit would 
support effective J7a(1) restoration. 



such, the degeneracy between pseudo-scalar and vector channels in the Wilson 
discretization [HEIEHI is not seen for staggered quarks. However, these results 
were not yet attempted to be extrapolated to the continuum limit and it is thus 
interesting that a recent paper jSJ reports that also in the latter discretization 
this degeneracy is reached in the continuum limit. 

Apart from screening masses, also the temperature dependence of spatial 
wave functions has been investigated [691 I77| . Also in this quantity, significant 
differences to zero temperature results could not be detected for temperatures 
T<T C . Above the critical temperature, the spatial wave functions become nar- 
rower with increasing temperature. Qualitatively, this observation is in accord 
|9U| with a rising spatial string tension |28l I89| . However, a quantitative com- 
parison is still lacking. 

Attempts to extract masses, rather than screening masses, from fits to tem- 
poral correlation functions have used extended operators of various kinds possi- 
bly combined with anisotropic lattices |fcjlUt>8ll75] . Both methods are means to 
counteract the short extent in the temperature direction. In the first case one 
hopes for a ground state contribution (if any) dominating already at small tem- 
poral distances. In the second case the isolation of the lowest lying state may 
be helped by the increased number of data points stabilizing more sophisticated 
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Figure 18: Temporal (left) pseudo-scalar correlation functions at temperatures 
below and above T c . Below T c the correlation function for one value of the quark 
mass is shown whereas above T c simulations have been performed in the chiral 
limit. The right hand figure shows the chiral extrapolation of the central value 
of the temporal pseudo-scalar correlator below T c . The bars denote value and 
error of the critical hopping parameter obtained from the chiral extrapolation 
of pseudo-scalar screening masses. 



fit ansatze which include more than a ground state contribution. However, 
one still has to rely on fit ansatze and, if possible, use the quality of the fit 
to distinguish between the various models. Moreover, modelling the operators 
introduces biases. The results have thus turned out to depend somewhat on 
the method. In general, they are at least qualitatively in agreement with an 
important two-quark cut contribution. In addition, the behavior of wave func- 
tions was interpreted to indicate meta-stable bound states |68|. In any case, the 
existence of genuine, narrow bound-states above T c could be excluded. 

6.3 Spectral functions from hadronic correlation functions 

The discussion in the previous section shows that we have evidence from numer- 
ical calculations that the spectral properties of hadrons change with tempera- 
ture. The drastic changes of hadronic correlation functions that occur when 
one crosses the transition temperature to the high temperature plasma phase 
are self-evident from the changes in the screening masses. The temperature 
dependence of correlation functions is particularly striking in the pseudo-scalar 
channel (see Fig. lT%|l . Below T c the central value, Gps(tT — 1/2, T), diverges in 
the chiral limit, which reflects that the pseudo-scalar remains massless fEa. l49|) . 
and above T c it stays finite and comes close to the free field limit. 

A much more subtle question is to quantify which thermal modifications 
of the spectral functions lead to the observed modification of the correlation 
functions. In order to decide on this we would like to have a more direct access to 
the spectral functions. It has been suggested recently [5T] to apply the Maximum 
Entropy Method (MEM), a well known statistical tool for the analysis of noisy 
data [55], also to the analysis of hadronic correlation functions. At least in 
principle, i.e. given sufficiently accurate numerical data on large lattices, this 
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approach allows to extract detailed information on hadronic spectral functions 
at zero as well as finite temperature. 

A lattice calculation of a hadronic correlation function on a lattice with tem- 
poral extent N T provides a set of data, {GH{Ti)jY i z\' 1 ^ D , where Nn denotes the 
number of measurements of a correlation function at a discrete set of Euclidean 
times TiT = i/N T , i = 1, ..,N T . The statistical analysis of this data set based 
on MEM aims at a determination of the most likely spectral function which 
describes the data given any prior knowledge on the correlation function or the 
spectral function. This information is taken into account in the so-called default 
model. In the case of hadronic correlation functions one usually tries to provide 
in this way information on the perturbatively known short distance behavior of 
Euclidean correlation functions. The first studies of spectral functions based on 
the MEM approach |)Tj and in particular the analysis of simple toy models at 
zero |91l 194] and finite temperature J23 have, indeed, been encouraging. The 
MEM approach has since then been used and further tested in various studies 
of hadronic correlation functions at zero and finite temperature both for light 
and heavy quarks [SI HD E3 EH EH1 11001 DUD GHH ■ We will discuss some of the 
results in the following. More details on the maximum entropy methods and its 
utilization for lattice QCD calculations are given in Ref. |93|. 

To illustrate the kind of thermal modifications of Euclidean correlation and 
spectral functions and the way this is analyzed using MEM let us briefly discuss 
some basic aspects of the correlation and spectral functions at zero and non-zero 
temperature. To be specific we consider the vector channel. The case of a freely 
propagating quark anti-quark pair is of relevance for the short distance part of 
the correlation functions as well as in the high (infinite) temperature limit. For 
massless quarks the corresponding spectral functions are given by Eq. 1511 

The corresponding correlation functions are shown in the left hand part of 
Fig. ^3 Note that the spectral functions for T = as well as T > are scale 
invariant, i.e. the rescaled spectral functions, ov/T 2 , depend on the rescaled 
energies, po/T, only. As a consequence the rescaled correlation functions Gy/T 3 
are also scale invariant and are functions of tT only. As a first test for the MEM 
analysis we may provide a set of data corresponding to the T > correlation 
function shown in Fig. ^5] (left) and use the T = spectral function as a default 
model. As can be seen in Fig. EH (right) the MEM analysis indeed can reproduce 
the thermal modification of the spectral functions already with a rather small 
set of equally spaced data points 95^ and, in turn, yields a perfect reconstruction 
of the correlation function. 

At zero temperature, the free spectral function (Fig. I19fl . only character- 
izes the short distance behavior of the correlation function or the large energy 
behavior of the spectral function, respectively. The most prominent modifica- 
tion, of course, results from the presence of the p-resonance and a more realistic 
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Figure 19: Correlation functions for a free quark anti-quark pair (left) and the 
corresponding spectral functions (right). The boxes indicate 32 equally spaced 
points on the free correlation function for T > 0. The curve through these 
points is obtained from the reconstructed spectral function. The upper curve in 
the left hand figure and the straight line in the right hand figure correspond to 
the input default model for the MEM analysis. 




spectral function is |104| . 

7T 



p2 r p m p 



+ 



p (pg — mp 1 ) 2 + r 2 m 
8ir\ 7r / 1 + exp((w -Po)/S) 



1 



(57) 



with T p = m p /(487rF 2 ) in the chiral (m* = 0) limit and m p = 770 MeV, 
F p = 142 MeV, lo = 1.3 GeV, 5 = 0.2 GeV and a s = 0.3. Due to the explicit 
appearance of hadronic scales this spectral function, of course, is no longer scale 
invariant in the sense introduced above. Assuming Eq. |23to hold also at non- 
zero temperature thus will also lead to vector correlation functions which are 
no longer scale invariant. In Fig. [2U] we show the scaled vector meson corre- 
lation functions, G p (tT)/T 3 , for temperatures T = 0.1, 0.2 and 1 GeV. The 
temperature dependence of the correlation function is apparent. At the high- 
est temperature G p (tT)/T 3 almost coincides with the corresponding rescaled 
free, zero temperature, correlation function. This shows that at these high tem- 
peratures the structure of the correlation function is dominated by the high 
energy part of the continuum contribution; the contribution of the p-resonance 
is suppressed. 

The scaling violations visible in Fig. 1201 are most prominent for tT ~ 0.5, 
which reflects that the hadronic scales, which violate the scaling of the spectral 
functions with temperature, show up at low energies. This is quite distinct 
from finite-cut off effects on the lattice that result in a distortion of the high 
energy part of the spectral functions and consequently lead to a modification of 
correlation functions at short distances. 

On finite lattices the temporal correlation functions built from free, massless 
quark propagators in the Wilson discretization jl()5| are given by a sum over 
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Figure 20: Vector meson correlation functions at temperatures below and above 
T c constructed from a phenomenological zero temperature vector meson spectral 
function. 



the quark momenta k, at zero spatial "meson" momentum 10 



G f ™ c ' lat (r,0)/T 3 = N c 



B) 2 cosh 2 (EN T /2) 



{a 1 ^ cosh[2^(r - N T /2)] + d 1 ^} 



(58) 



Here, B is given by B = 2 J^. sin 2 (fci/2) and E denotes the quark energy E — 
|^/a74+ y/l Ta/il with a = (A + B 2 )/(l + B) where A = ^ 



2 In 



sin 2 kj 



The coefficients and are collected in Table Q] for some channels H. 
Note that the coefficient d}^ approaches 1 in the continuum limit. Moreover, in 
the zero temperature limit, the r independent piece of the correlation function 
G^ ee ' lat (r, 0) vanishes proportional to T 3 . In Fig. [2] we show the free vector 
correlation functions and the corresponding spectral functions |ll)3j on lattices 
with temporal extent N T — 16, 24 and 32. As can be seen the influence of 
finite cut-off effects is restricted to the short distance behavior of the correlation 
functions, which results from the strict momentum cut-off on a finite lattice. 
We also note the pronounced peak in the lattice spectral function, which results 
from the distortion of the lattice dispersion relation for free quarks at large 
momentum and, in particular, close to the corners of the first Brillouin zone. In 
the interacting case peaks at apo ~ 2 thus have been attributed to the influence 
of the heavy Wilson doublers [51] . 



10 Similar formulae are available for staggered quarks 1781 . 
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Figure 21: Free pseudo-scalar correlation functions (left) and spectral functions 
(right) calculated on lattices with temporal extent N T using the Wilson fermion 
formulation. 



6.4 Spectral analysis of thermal correlation functions 

In Section 6.3 we have discussed properties of the pseudo-scalar correlation 
functions below and above T c which clearly suggest that the spectral properties 
in this channel drastically change when going from the low temperature hadronic 
phase to the high temperature plasma phase. After this drastic change the 
rescaled pseudo-scalar correlation function above T c , however, seems to be only 
weakly temperature dependent. This property is even more pronounced in the 
vector channel. In Fig. [221 we show the vector correlation function for a set 
of temperatures ranging from 0.4 T c up to 3 T c calculated on lattices of size 
up to 64 3 x 16. As can be seen the rescaled correlation functions are within 
statistical errors temperature independent. In view of the discussion presented 
in the previous section this is indeed remarkable; the scaling of the correlation 
function is inconsistent with a temperature independent spectral function and, 
in fact, suggests that also the parameters determining the spectral function (or 
at least its dominant part) scale with temperature. As at high temperature 
the vector spectral function is dominated by the continuum contribution the 
most obvious expectation is that the parameters controlling the threshold in 
the continuum contribution, uig and 6, increase with temperature. 

The existence of a temperature dependent energy cut-off in the vector as well 
as the pseudo-scalar spectral functions above T c is also obtained in the MEM 
analysis of the corresponding correlation functions. This is shown in Fig. [23 In 
both cases the reconstructed spectral functions suggest that the contributions 
from the pion and rho states disappear above T c . Instead, a broad peak becomes 
visible which shifts to larger energies with increasing temperature. 

6.5 Vector meson spectral function and thermal dilepton 
rates 

The vector spectral function shown in Fig. [23 is directly related to the thermal 
cross section for the production of dilepton pairs at vanishing momentum, 
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Figure 22: Vector meson correlation functions at various temperatures below 
and above T c . For better visibility data points have slightly been shifted hori- 
zontally relative to the data set at T — 0.4T C . 
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This thermal dilepton rate is shown in Fig. The "resonance" like en- 
hancement seen in Fig. [531 results here in the enhancement of the dilepton rate 
over the perturbative tree level (Born) rate, 
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dp a d 3 p 



\p=o 
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(60) 



for energies Pq/T £ [4, 8]. Obvious discrepancies with hard thermal loop calcu- 
lations |l()f>) show up at smaller energies where the lattice results for the spec- 
tral functions rapidly drop while the HTL result diverges in the infrared limit. 
In fact, this discrepancy exists already on the level of the spectral functions 
themselves. The lattice spectral functions vanish (rapidly) in the limit po — * 
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Figure 23: Pseudo-scalar (left) and vector (right) spectral functions obtained 
from quenched QCD calculations on lattices up to size 64 3 x 16 using improved 
Wilson fermions. 
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Figure 24: Thermal dilepton rate 

whereas in the vector channel the HTL-spectral function diverges As we 
have pointed out in the discussion of the p-meson spectral function at T = 0, 
with increasing temperature it becomes more and more difficult to resolve the 
low energy part of spectral functions. The suppression of the thermal dilepton 
rate at low energies deduced from the structure of oy thus will require further 
detailed investigations on large lattices. However, even without invoking the 
MEM analysis to determine ay (po , T) one can deduce a minimal constraint on 
the low energy behavior of the spectral function from the fact that the corre- 
lation function Gy(r, 0) stays finite at high temperature. This demands that 
a v{P0iT) ~ Pq with a > 0. The vector spectral function thus has to vanish in 
the limit po — > 0. In fact, in order to get non-vanishing transport coefficients in 
the QGP the spectral function should be proportional to po in this limit jl()7| . 
This may be included in a MEM analysis as an additional constraint |l()8j . In 
any case, getting the low energy behavior of the vector spectral functions under 
control still is a challenge for future studies. 

6.6 Heavy quark spectral functions and charmonium sup- 
pression 

As the properties of mesons constructed from light quarks are closely related to 
chiral properties of QCD it is expected that these states dissolve in the QGP. The 
situation, however, is different for heavy quark bound states, which generally 
are expected to be controlled by properties of the heavy quark potential at 
moderate distances (R<0.5fm). Moreover, this distance scale will depend on 
the quark flavor under consideration (charmonium, bottonium). Heavy quark 
bound states thus are sensitive to screening of the heavy quark potential at 
high temperature and therefore probe the deconfining aspect of the QCD phase 
transition [SHI - Whether a heavy quark bound state survives the QCD phase 




Thermodynamics and in-medium hadron properties from lattice QCD 



41 




4 8 12 16 



Figure 25: Heavy quark spectral functions and screening masses (vertical lines) 
in the vector channel calculated in quenched QCD on lattices of size 40 4 (T ~ 
0.9T C ), 64 3 x 24 (T ~ 1.5T C ) and 48 3 x 12 (T ~ 3T C ) at (3 = 7.19. 



transition or not strongly depends on the efficient screening of the interaction 
among quarks and anti-quarks in the QGP. Potential model calculations indicate 
that some cc bound states, e.g. the J/tp and in particular the bb bound states 
could survive at temperatures close to T c while radially excited states like the 
Xc most likely get dissolved at T c |109| . 

A more direct analysis of the fate of heavy quark bound states again can 
come from the analysis of thermal meson correlation functions constructed with 
heavy quarks |l()()l ll()2j . These first, exploratory studies show that pseudo- 
scalar (jjc) and vector (J/tp) correlation functions show almost no temperature 
dependence across the phase transition and can survive in the high temperature 
phase at least for temperatures T<l.bT c . In Fig. [211 we show results for the 
vector meson spectral function at temperatures below and above T c . They are 
compared with results for screening masses which are indicated by horizontal 
lines in this figure. As can be seen, the pole masses, identified as the location of 
the first peak in the spectral functions, coincide with the screening masses up to 
T = 1.5T C . Whether the apparent broadening shown in the figure and reported 
in Ref. jl()2| for pseudo-scalar and vector meson resonances is premature remains 
to be seen. At 3T C on the other hand there is no evidence for a pole in the 
spectral function anymore and the threshold in the spectral function is unrelated 
to the value of the screening mass which itself moves away from its value at 
the lower temperatures. Thus there is evidence for modifications of spectral 
properties at high temperature and there are also first hints that the vector 
meson resonance is dissolved at least at T ~ 3T C [TTHj . 
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7 Summary 

The survey of recent lattice calculations presented here shows that there has 
been significant progress in our understanding of the thermodynamics of strongly 
interacting matter since this field has been reviewed for the last time in this se- 
ries of publications |112j . We now have detailed quantitative information on the 
QCD equation of state at high temperature and the transition temperature to 
the plasma phase. Moreover, there exist now first calculations for non- vanishing 
chemical potential that allow us to judge the influence of a non-zero baryon num- 
ber density on the transition temperature and the equation of state. Although 
these calculations are at present restricted to the region of high temperature 
and small values of the chemical potential they give first insight into the QCD 
phase diagram in the T — /i plane and even gave first evidence for the presence 
of a critical point in this plane which has been anticipated in phcnomenological 
models. Nonetheless, lattice calculations still are performed with quark masses 
larger than those realized in nature and we have to proceed to smaller quark 
masses to provide reliable input to the analysis of current and future heavy ion 
experiments. 

One of the outstanding problems in QCD thermodynamics is to gain a quan- 
titative understanding of thermal modifications of hadronic properties at high 
temperature. Although it is generally expected that hadronic bound states can- 
not survive in the high temperature plasma phase we have little quantitative 
information about their dissolution in a hot thermal medium. Lattice calcu- 
lations have provided many indirect information on thermal modifications of 
hadron properties, e.g. through the calculation of hadronic screening lengths 
and susceptibilities. However, it was only recently that we gained first insight 
into thermal modifications of hadron masses through the application of statis- 
tical tools like the maximum entropy methods. Hopefully, this will lead in the 
future to a direct verification of the disappearance of light and heavy quark 
bound states, the structure of quasi-particle excitations in the plasma phase 
as well as experimentally and phcnomcnologically important observables like 
dilepton rates or even transport coefficients. 

Progress in lattice studies of the QCD phase diagram as well as the calcu- 
lation of hadron properties have significantly been influenced through new con- 
ceptual ideas and the application of new calculational schemes. Nonetheless, 
progress in lattice calculations has always been closely related to the progress 
made in computer technology. Very soon a new generation of computers with a 
sustained speed in the teraflops range will become widely available to the lattice 
community. It is to be expected that this will allow us to perform much more 
realistic numerical calculations with parameters that reflect the quark mass 
spectrum realized in nature. 
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A Appendix 

In this appendix we collect some formulae hopefully useful in the context of our 
discussion of hadronic correlation functions. Full account of the imaginary time 
formalism (at finite temperature) is given in textbooks |71 1591 IfiH] . 

For the generic case of a spinless Bose field the Feynman propagator in real, 
Minkowski time t is given by 

D(t,x;t',x l ) = {T{$(t,x)^(t',x')}) . (61) 

Here, the expectation value is taken at finite temperature, 

(T4(t,x)ft(t',tf )> = l-tTTj>(t,x)ft(t',x')e- ii/T , (62) 

with H being the Hamiltonian and Z the partition function. The Feynman 
propagator can be written as the sum of two terms 

D(t, x; t', x') = Q(t - t')D> (t, x; t', x") + 9(f' - t)D<(t, x; t', x') (63) 

where 

D > {t,x;t',x") = {cj>(t,x)$(t',x*)) = D < {t',x';t,x) . (64) 
The cyclicity of the trace in Eq. leads to the KMS condition, 

D>(t;t') = D<{t + i/T;t') . (65) 

In fact, the functions D, D > and D < only depend on coordinate differences 
because of translational invariance. Thus, their Fourier transformation is given 
as, for instance, 

D>(t-t',x-x>)= f £L [ + ™ ^e- i ^-*'^-^D>(p ,?) . (66) 
J (27r) d J^oo 2vr 

The spectral density cr(po,p) is defined as 

a(po,p)=2w[D>(p ,p)-D<(p ,p)] . (67) 

This also shows that o~(po,p) is an odd function of the energy variable pq. Note 
that the spectral density not only depends on the Lorentz invariant p 2 as at 
zero temperature but also on the scalar product np with = (1,0) which 
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arises from the presence of the heat bath. By means of Eas.l64landl66lit is also 
obtained as 

a(p ,p) = 2tt / d 3 x / dt e^ ot ~^ ( [$(t, x),^(0, 0)] ) . (68) 



Eq. E3 together with the KMS relation can be inverted to give 

^^= l- e - P0 /T 2n ■ (69) 

In imaginary time, t = —ir with r real, one can introduce a "time" ordering 
in t and define 

A(t,x;t',x , ) = (T{^(t,x^(t',x')}) . (70) 

The periodicity condition now follows from Eq. 1651 

A(t,£) = A(t,£;0,0) = A(t-1/T,£;0,0) . (71) 

As a consequence the Fourier spectrum consists of discrete Matsubara frequen- 
cies, 

w„ = 2wnT,n = 0,±1,±2, ... , (72) 
at which the Fourier transform is defined 

f f 1/T 

A(ioj n ,p) = / d 3 x / dTe^ T ~ %px A(r,f;0,0) . (73) 



The inverse Fourier transform then includes a sum over Matsubara frequencies 
A(T,x)=TY,J-^e- iuJnT+ipS A{iu n ,® . (74) 

Due to the identification 

A(r) = D>(-i T ) (75) 
in the r interval [0, 1/T] one then obtains from Eas. 1731 and 1691 

A(^,p)= r dpQ ^L . (76) 

J-oc PO - l^n 

Under certain assumptions, this expression can be analytically continued to 
complex frequencies z. In fact, for values z = fco ± it with kg real, one obtains 
the retarded and advanced propagators, 

D R {t) = e(*)(0(t)^(o)) , 

D A (t) = -Q(-t)(i(t)ft(p)) , (77) 

as 

D R/A (k )=TA(k ±ie) . (78) 
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On the lattice one usually studies the temporal correlator at fixed momentum 
p. By doing the sum over Matsubara frequencies this can easily be obtained 
from Eqs. [To] and [7i>| as 



OrlrM - j (Fx A(r,x)e~ 



"dp„ <r(p„,p) e-"'[e(T)+n B (p )] (79) 

J — oo 

with the Bose distribution 

n B (po) = [exp(p /T) - l]- 1 . (80) 
Without loss of generality we take r positive, re[0, l/T], to arrive at Eas. 1361 

Instead of the temporal correlator, most (lattice) analyses have been concen- 
trating on the spatial correlation functions. These depend of course on the same 
spectral density but are different Fourier transforms of it in r and x± — (x, y) 

f 1/T f 

G(ioj n ,p±,z) = I dr d 2 x± exp[zo; ra T — ip±x±] A(r, x) 



o J 

+ °° dPz lp , z f + °° , 2po (j(p ,p±,Pz) , Q1 , 
— e"' / dp — — - (81) 

Projecting onto vanishing transverse momentum and vanishing Matsubara fre- 
quency then gives Eq. |2H1 

For a free stable boson of mass M the spectral density has been given by 
Eq. 1421 Setting the matrix elements to 1 , in Euclidean momentum space the 
propagator is thus obtained as 

A(iu n ,p) = 2 l (82) 

w ; + 

Computing the correlator in imaginary time gives 



Gt(t,p) 



2uj p 1 - e -"vl T 

^K(u p ,t) (83) 



where 



lu p = V? 2 + M 2 (84) 



In the limit of vanishing momentum the correlator thus decreases with the 
mass modulo periodicity, Eq. 1521 On the other hand, the spatial correlator is 
computed as 

G(uj n ,p ± ,z) = ^—e-^ z (85) 
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where w sc is given as 

u sc = yJu*+pl+AP (86) 

The opposite limit is reached for two freely propagating quarks contributing 
to the spectral density. Here the starting point is the momentum space rep- 
resentation of the meson correlation functions. To leading order perturbation 
theory one has to evaluate the self-energy diagram shown in Fig. 114b , in which 
the internal quark lines represent a bare quark propagator Sp{iijj ni p) [S3J. The 
quark propagator can be expressed in terms of its spectral function 11 <7f roc (p , p) , 

CTfroo^cP) = e(p ) {i> + m) S(pI - p 2 - m 2 ) (87) 
and is conveniently written as 

S F {iu} n ,p) = / dpo : (88) 

J-oo Po - iu n 

Contrary to Eq. [721 the fermion Matsubara frequencies are odd integer multiples 
of ttT, 

w n = (2n + 1)ttT (89) 

A free thermal meson spectral function is then obtained from the imaginary 
part of the correlation function in momentum space, confer Eq. 1761 

A H (iuj m ,p) = N C T*^2 J (2^)3 tr ^HSF(iun, k)F' H Sp(iuj m - iu> n ,p- k) 

(90) 

where uj m denotes a boson frequency and lu ti a fermion one. The matrices Th 
ensure projection on the spin and parity quantum numbers of a mesonic state 
H. 

The above analysis for the free thermal quark-antiquark correlators can be 
extended to the leading order HTL approximation. In this way important 
medium effects of the quark-gluon plasma such as effective quark masses and 
Landau damping are taken into account. The HTL-resummed fermion prop- 
agator is obtained from Eq. 1881 by replacing the free spectral function <7f ree 
with the HTL-resummed spectral function which for massless quarks is given 

by nnnm 

o-htl(po,p) = ^<?+{po,p){jo -Pi) + ^-{po,p)ilo +Pl) (91) 
with p = p/p, p = \p\-, and 

2 2 

v±{Po,p) = -|^2— i s (po - w ±) + s (Po + + f3±(p ,p)O(p 2 -p ) 



n We use the convention as in Ref. 1591 where 74 = 170 and po — * —ipi 
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P±(po,p) = -^(PTPoK PiPTPo) + m 2 T ( 1 ± 



PTPo 
2p 



In 



P -Po 



P + Po 



) 



r 7T 2 pTPo] 2 \ 
+ l2 mT —\ ) 



(92) 



Here u>± (p) denote the two dispersion relations of quarks in a thermal medium 
|1U61 llll| and nix — g(T)T/y/6 is the thermal quark mass. In addition to 
the appearance of two branches in the thermal quark dispersion relation the 
HTL-resummed fermion propagator also receives a cut-contribution below the 
light-cone (pq < p 2 ), which results from interactions of the valence quarks with 
gluons in the thermal medium (Landau damping). Furthermore, an explicit 
temperature dependence only enters through mx(T). It also should be noted 
that the HTL resummed quark propagator is chirally symmetric in spite of the 
appearance of an effective quark mass 



